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ABSTRACT 

O 

^ ; We present a catalog of 59 z ~ 4.5 Lya-emitting galaxies spectroscopically 

confirmed in a campaign of Keck/DEIMOS follow-up observations to candidates 
selected in the Large Area Lyman Alpha (LALA) narrow-band imaging survey. 
We targeted 97 candidates for spectroscopic follow-up; by accounting for the 
variety of conditions under which we performed spectroscopy, we estimate a se- 
lection reliability of ~ 76%. Together with our previous sample of Keck/LRIS 
Oh 1 confirmations, the 59 sources confirmed herein bring the total catalog to 73 spec- 

troscopically confirmed z ~ 4.5 Lya-emitting galaxies in the ~ 0.7 degrees 
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covered by the LALA imaging. As with the Keck/LRIS sample, we find that 
a non- negligible fraction of the confirmed Lya lines have rest-frame equivalent 
widths (W^ Bt ) which exceed the maximum predicted for normal stellar popula- 
tions: 17% - 31% (93% confidence) of the detected galaxies show W{ cst > 190 A, 
and 12% - 27% (90% confidence) show W{ cst > 240 A. We construct a luminosity 
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function of z ~ 4.5 Lyct emission lines for comparison to Lya luminosity func- 
tions spanning 3.1 < z < 6.6. We find no significant evidence for Lya luminosity 
function evolution from z ~ 3 to z ~ 6. This result supports the conclusion 
that the intergalactic medium remains largely reionized from the local universe 
out to z Rj 6.5. It is somewhat at odds with the pronounced drop in the cos- 
mic star formation rate density recently measured between z ~ 3 and z ~ 6 in 
continuum-selected Lyman-break galaxies, and therefore potentially sheds light 
on the relationship between the two populations. 

Subject headings: cosmology: observations — early universe — galaxies: evolu- 
tion — galaxies: formation — galaxies: high-redshift 



1. Introduction 



Observational cosmology has recently witnessed a tremendous increase in proficiency 
in the identification of galaxies at the earliest cosmic epochs. Thanks in large part to the 
availability of large-format mosaic CCDs well-suited for wide-field imaging and spectroscopic 
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')) to the assembly of statistically robust samples spanning the earli- 



est accessible redshifts. Robust samples of this kind are necessary for understanding the 
systematics of selection criteria, and of the spatial distribution of the galaxies themselves. 
Deficiencies in such understanding are the main source of uncertainty in inferred luminosity 
functions and universal star formation rates, which in turn are the keys to understanding 
the cosmic history of star formati on, galaxy assembly and evolution, and eve n the early 



ionization history of the IGM (e.g., Malhotra fc Rhoadd 12004 ; IStern et all 120051 ) 
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2004al B I Yan Sz Windhorstl I2OO4J ) . These two techniques are complementary; Lya searches 
at typical sensitivities can identify galaxies with UV-continua too faint to be detected by 
the Lyman break method, but such surveys only select that fraction of galaxies with strong 
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line emission. 



The Large Area Lyman Alpha (LALA) survey (jRhoads et al.ll2000l ) has recently iden- 
tified in deep narrow-band imagin g a large sample of Lya-emitting galaxies at redshift s 
z « 4.5 jMalhotra & Rhoadsll2002h . z re 5.7 JRhoads fe Malhotrall200ll ; iRhoads et al.ll2003h . 



and z ~ 6.5 (jRhoads et alihopj ). In bawson et all I2004J (Paper I), we reported on the 



spectroscopic confirm ation with the W. M. Keck Observatory's Low Resolution Imaging 
Spectrometer (LRIS; Oke et all 1 19951 ) of 17 Lya-emitting galaxies selected in the LALA 
z ~ 4.5 survey. The resulting sample of confirmed Lya emission lines showed large equiv- 
alent widths (median W{ CBt ~ 80 A) but narrow velocity widths (FWHM Av < 500 km 
s -1 ), indicating that the Lya emission in these sources derives from star formation, not from 
AGN activity. Models of star formation in the early universe predict that a small fraction of 
Lya- emitting galaxies at z ~ 4.5 may be nascent, metal-free objects (e.g.. lScannapieco et al. 
20031 ) . and indeed we found with 90% confidence that 3 to 5 of the confirmed sources exceed 
the maximum Lya equivalent width predicted for normal stellar populations. However, we 
did not detect the He II A1640 emission expected to be characteristic of primordial star for- 
mation. Specifically, the He II A1640 flux in a composite of the 11 highest resolution spectra 
in the Keck/LRIS sample was formally consistent with zero, with a 2a (3a) upper limit of 
13% (20%) of the flux in the Lya line. In other words, though these galaxies may be young, 
they show no evidence of being truly primitive, Population III objects. 

We have recently more than quadrupled our catalog of spectroscopically-confirmed Lya- 
emitting galaxies at z ~ 4.5 with a spectroscopic campaign u sing the W. M. K eck Obser- 
vatory's DEep IMaging Multi-object Spectrograp (DEIMOS; iFaber et al.ll2003l ). targeting 
candidates selected from the LALA survey. Together with the detections presented in Pa- 
per I (and accounting for minor overlap in the samples), the 59 Lya-emitters confirmed 
with Keck/DEIMOS bring the total catalog to 73 spectroscopically confirmed z rs 4.5 Lya- 
emitting galaxies in the ~ 0.7 degrees 2 imaged by LALA. In this paper, we utilize these 
additional confirmations to update the results of Paper I, and to construct a luminosity 
function of z ~ 4.5 Lya emission lines for comparison to Lya luminosity functions spanning 
3.1 < z < 6.6. We describe our imaging and spectroscopic observations in § [2], and we sum- 
marize the results of the spectroscopic campaign in §[31 In §Hl we investigate the distribution 
of the Lya lines in equivalent width, we construct Lya luminosity functions for our sample 
and for several extant samples, and we discuss the implications of the luminosity functions 
for the relationship between Lya-emitters and Lyman-break galaxies (LBGs), and for the 
history of reionization. Throughout thi s paper we adopt a A-cosmology with Q M = 0.3 and 
= 0.7, and H = 70 km s _1 Mpc -1 (ISpergel et al.l 120031 ) . At z = 4.5, such a universe is 



1.3 Gyr old, the lookback time is 90.2% of the total age of the Universe, and an angular size 
of l'.'O corresponds to 6.61 comoving kpc. 
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Observations 



2.1. Narrow-Band and Broad-Band Imaging 



The LALA survey concentrates on two primary fields, "Bootes" (14:25:57 +35:32; 
J2000.0) and "Cetus" (02:05:20 -04:55; J2000.0). Each field is 36 x 36 arcminutes in size, 
corresponding to a single field of the 8192 x 8192 pixel Mosaic CCD cameras on the 4m 
Mayall Telescope at Kitt Peak National Observatory and on the 4m Blanco Telescope at 
Cerro Tololo Inter- American Observatory. The z « 4.5 search uses five overlapping narrow- 
band filters each with full width at half maximum (FWHM) 80 A (Figure [1]). The central 
wavelengths are 6559, 6611, 6650, 6692, and 6730 A, giving a total redshift coverage of 
4.37 < z < 4.57 and a survey volume of 7.4 x 10 5 comoving Mpc 3 per field. In roughly 
6 hours per filter per field, we achieve 5a line detections in 2'.'3 apertures of ~ 2 x 10~ 17 
erg cm~ 2 s _1 . 

The primary LAL A survey fields were chosen to lie within the NOAO Deep Wide- 
Field Survey (NDWFS: IJannuzi fc Deylll999l ). Thus, deep NDWFS b road-band images are 
available in a custom B w filter (A = 4135 A, FWHM = 1278 A; l.lannuzi fc Devlll999 



Jannuzi et al., in preparation) and in the Harris set Kron-Cousins R and J, as well as J, 
H, K, and K s . The LALA Bootes field benefits fro m additional deep V and SDSS z' filter 
imaging. The imaging data reduction is desc r ibed i n lRhoads et al . (2000) , and the c andidate 
selection is described in lRhoads &: Malhotral ( J200lh and lMalhotra fc Rhoadd (120021 ). Briefly, 
candidates are selected based on a 5a detection in a narrow-band filter, the flux density of 
which must be twice the i?-band flux density, and must exceed the i?-band flux density at 
the 4<t confidence level. To guard against foreground interlopers, we require an observed 
equivalent width W^ hB > 80 A, and a non-detection in the B w band (at the < 2a level). 



2.2. Spectroscopic Observations 

Between 2003 March and 2004 May we obtained spectroscopy of 97 z ~ 4. 5 candidate 



Lyq- emitters with the DEep IMaging Multi-object Spectrograph (DEIMOS; iFaber et al. 



20031 ). a second-generation camera on the Keck II telescope with high multiplexing capa- 
bilities and improved red sensitivity. Each slitmask included approximately 15 candidate 
Lya-emitters (mixed in with roughly 50 other spectroscopic targets) and was observed for 
1.5 to 2.0 hrs in 0.5 hr increments. Six slitmasks targeting a total of 80 candidates were 
observed in the Bootes field; the airmass in these observations never exceeded 1.5. One slit- 
mask targeting 17 candidates was observed in the Cetus field; the airmass for this slitmask 
was constrained to less than 1.8. The seeing in all observations ranged from 0'.'5 to l'.'O. We 
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estimated the seeing by examining the alignment stars observed during the direct-imaging 
phase of setting up the slitmask; that is, the stars were imaged through the BAL12 (clear) 
filter with the grating angle set such that zero-order light fell on the detector. 

All observations employed l'.'O wide slitlets and the 600ZD grating (Abi aze = 7500 A; 
0.65 A pixel -1 dispersion; AAfwhm ~ 4.5 A ps 200 km s _1 )0. The wavelength range covered 
by a typical slitlet was roughly 5000 A < A < 10000 A. The precise wavelength coverage 
depended somewhat on the location of the slitlet on the slitmask, but never exceeded 4390 
A at the low extremum or 1.1 /mi at the high extremum. No order-blocking filter was used; 
since the targets were primarily selected to have red colors, second order light should not 
be of concern. Most nights suffered from s ome cirrus; relative flux ca libration was achieved 
from observations of standard stars from iMassey fc Gronwalll (119901 ) observed during the 
same observing run. It should also be noted that the position angle of an observation was 
set by the desire to maximize the number of targets on a given slitmask, so observations 
were generally not made at the parallactic angle. 

We processed the two-dimensional data using the DEEP2 DEIMOS pipeline^. We per- 
formed small (0 / /5) dithers between exposures on our initial observing run; to reduce these 
data, we supplemented the DEEP2 DEIMOS pip eline with ad ditional home-grown routines. 
We extraeted^eetra with the IRAF@ package (iTodyl Il993l ) using the optimal extraction 
algorithm (jHornd ll986). following standard slit spectroscopy procedures. 



Prior experience with faint-object spectroscopy dictates that a small but significant er- 
ror in the measured flux of faint continua may be introduced by sky subtraction during the 
processing of the two-dimensional spectra. We investigated this possibility with ~ 10 addi- 
tional non-overlapping extractions in source-free regions in each two-dimensional spectrum, 
parallel to and along the same trace as the extraction of the neighboring Lya-emitting galaxy. 
We then fitted these blank-sky spectra over the same region that we fitted for the contin- 
uum redward and blueward of the emission line in the object extraction. For 15 sources, 
these fits yielded a tiny residual signal, which we interpreted as a systematic error in the 



1 We measured the instrumental resolution by autocorrelating one-dimensional extracted spectra of night- 
sky emission lines. The autocorrelation results in an effective average line profile with a high signal-to-noise 
ratio, which we fit with a Gaussian to obtain the FWHM. We performed this test on ~ 50 night-sky spectra 
with the result AAfwhm = 4.47 ± 0.03 A. The quoted uncertainty is the error in the mean, and does not 
include possible systematic effects due to blended night sky lines. 

2 See http : / / astron . berkeley . edu/~cooper/deep/spec2d/ 

3 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Associa- 
tion of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 
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t wo- dimensional sky subtraction and applied as a correction to quantities derived from the 
object spectra. The typical correction was ~ 0.04 ± 0.02 /xJy, but the correction reached as 
high as ~ 0.1 ± 0.07 /iJy in three cases. Sky-subtraction residuals of this kind generally re- 
sulted when a small spectroscopic slit contained a bright serendipitous detection in addition 
to the target, the combination of which made it difficult to fit the sky background. 

3. Spectroscopic Results 

Out of 97 spectroscopic candidates, we achieved 73 detections, 59 of which constitute 
Lya confirmations according to the criteria outlined below. A histogram of these confir- 
mations appears in Figure [TJ and a set of sample spectra are shown in Figure [2j The 
spectroscopic properties of the Lya confirmations are summarized in Table [D One detected 
galaxy lacks an emission line but shows a large spectral discontinuity identified as the onset 
of foreground Lya-forest absorption at z = 4.4620. Three of the detections are identifi- 
ably low-redshift interlopers (two resolved [O II] A3727 doublets at z ~ 0.8; one complex of 
[O III] AA4959,5007 and H/5 at z ~ 0.3) which survived the candidate selection thanks to 
their unusually high equivalent widths (e.g., W^ s > 2000 A). In ten cases, we see a possible 
low signal-to-noise ratio (< 1) emission line located at the correct location in both wave- 
length space and physical position to be associated with the narrow-band selected target. 
However, even if these "detections" are real, they cannot be reliably identified as either Lya 
or as low-redshift interlopers. If these 10 cases were in fact low signal-to-noise ratio detec- 
tions of Lya emission, then the "success rate" of the Keck/DEIMOS campaign would be 
72%, identical to that of the Keck/LRIS sample described in Paper I (but also subject to 
all the caveats listed therein). We do not include unconfirmed sources in any of the ensuing 
discussion. 

The remaining 24 targets were classified as nondetections. Five of these slitlets suffered 
from some kind of instrument or reduction issue, e.g., the target was dithered off the slitlet 
and so did not reproduce across the individual integrations, or irregularities in the machining 
of the slitmask resulted in defects in the data processing. Of the final 19 nondetections, 13 
targets were observed under adverse conditions (e.g., variable cloud coverage, and/or poor 
seeing) for which the general spectroscopic yield was low. Our failure to confirm these targets 
as z ~ 4.5 Lya-emitters should not be taken to bear on the efficacy of candidate selection. 



A sufficiently bright Lyman break galaxy can be selected as a narrowband excess object when the narrow 
filter lies redward of the Lya forest, so that neutral hydrogen absorption significantly reduces the broad band 
flux without affecting the narrow band flux. 
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Six nondetections were observed under photometric conditions with subarcsecond seeing 
for which the spectroscopic yield was otherwise high. However, subsequent inspection of the 
imaging revealed that five of these targets were suboptimal candidates for one of a variety 
of reasons: two candidates sit on weak satellite trail residuals; one candidate appears in 
an initial epoch of imaging but not in subsequent epochs, suggesting that it is a variable 
source or a spurious detection; two candidates are marginal or irregular detections in the 
imaging. This leaves just one otherwise viable candidate Lya-emitter that was not confirmed 
in spectroscopy, even though the conditions for spectroscopy were favorable. Since this source 
(J1424398+353801) was a single-band detection in the narrow-band imaging, it is possible 
that it represents a spurious false-positive and not a genuine candidate. Given the large 
number (10 6 5 ) of independent resolution elements in the images, we expect about one false- 
positive at the 5<7-level per LALA field per narrow-band filter, and t his number could b e 



larger if the noise properties of the image are not precisely Gaussian (see lRhoads et al.ll2003l ). 



To estimate the reliability of our candidate selection, we consider only the foregoing six 
non-detections, the three low-redshift interlopers, and the 10 low signal-to-noise ratio detec- 
tions as legitimate non-confirmations. This admittedly rough scheme suggests a rate of 59 
detections out of 78 viable candidates observed spectroscopically under workable conditions, 
for a final selection reliability of ~ 76%. The rate of spectroscopic confirmation is plotted 
as a function of narrow-band flux in Figure [31 



3.1. Spectroscopic Sensitivity to Lya Emission 

We now assess our spectroscopic sensitivity to line emission. Each one-di mensional 



spectrum was created with a variance-weighted optimal extraction (lHorndll986l ) from the 
two-dimensional data. For each object, we therefore have both the flux and the flux variance 
as a function of wavelength. We used the variance spectrum to estimate the uncertainties in 
quantities derived from the object spectra, e.g., we sum the variance spectrum in quadrature 
over the wavelength range covered by the observed line profile to estimate the uncertainty 
in the measured line flux. However, the variance may also be used to measure the noise over 
wavelength ranges corresponding to any Lya line we might have detected given the redshift 
range permitted by our narrow-band imaging, roughly 4.37 < z < 4.57. Accordingly, for 
each object we ranged over redshift and calculated the smallest emission line flux detectable: 



Aa(«) 

F Um (z) = n sie 5 disp I ^2 °\ 
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where n B i g is the minimum signal-to-noise ratio necessary for a detection (here taken to be 
3), <5 disp is the grating dispersion (0.63 A pix -1 for the Keck/DEIMOS 600ZD grating), and 
&\ is the flux error in each pixel determined during the variance-weighted one-dimensional 
extraction, in units of fx- The limits Ai and A2 are defined by 

A 1 (2) = (l + ^)(1216-AA/2), 

X 2 (z) = (l + z)(1216 + AA/2) , 1 ' 

where A A is the fiducial rest-frame full width of the emission line (here taken to be 3 A). 

We assembled the F lim (z) for each object into a grid and then ranked the F lim at each 
redshift, resulting in the cumulative distribution of sensitivity to Lya emission line flux 
shown in Figure HI The distribution may be interpreted as giving the probability that a 
putative Lya emission line of a given flux and a given redshift would have been detected in 
our spectroscopic campaign. Since we cover a comparatively small redshift range centered 
essentially at the peak of the detector throughput, the sensitivity distribution is dominated 
entirely by night-sky emission lines rather than by instrumental effects. And since the 
original narrow-band survey was designed to probe relatively noise-free windows in night- 
sky emission, the spectroscopic sensitivity is fairly flat over the redshift range of interest. 
In sum, the implied depth of our spectroscopic survey is 50% complete to /(Lya) ~ 3 x 
10~ 18 erg cm -2 s^ 1 , approximately 7 times deeper than the narrow-band imaging. Note 
that because we derived this sensitivity function from the sample of spectra themselves, it 
depends entirely on the details governing the manner in which these spectra were obtained 
and processed, and is therefore valid for this survey only. 



3.2. Redshift Identification 



Of course, given the detection of an emission line, the identification of that line as 
high-redshift Lya can remain problematic. Thor ough treatments of the pitfalls of one-line 



redshift identificatio ns are given elsewhere (e.g., IStern fc Spinradl Il999l : IStern et al.l l2000l ; 
Dawson et al.ll200ll ). In surveys of the present kind, the primary threat to the proper in- 



terpretation of a solo emission line is the potential for low-redshift, high-equivalent width 
[O II] A3727 to survive candidate selection, and then to be misidentified as high-redshift Lya; 
in later spectroscopy. However, at z = 0.8 (the redshift of an [O II] A3727 line mistaken for 
Lya at z = 4.5), the redshifted separation between the individual lines of the [O II] A3727 
doublet (rest wavelengths 3726 A and 3729 A, respectively) is 5.4 A. The doublet is therefore 
just resolved in our spectroscopy and serves to uniquely flag [O II] A3727 interlopers (Fig- 
ure [5]); this is an improvement afforded by Keck/DEIMOS over the spectroscopy presented 
in Paper I. Less frequently, high-equivalent width [O III] A5007 survives as an interloper in 
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our candidate selection. However, [O III] A5007 can typically be identified by neighboring 
[O III] A4959 at one-third its strength, or by neighboring B./3. 

Beyond merely eliminating plausible low-redshift interlopers, we may identify Lya emis- 
sion by its characteristically asymmetric morphology, or by the presence of a continuum break 
if the continuum is sufficiently well-detected. Each of our confirmed Lya detections demon- 
strates the asymmetric emission line profile characteristic of the line, where neutral hydrogen 
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present a scatter plot of the flux-based asymmetry statistic: 

A A r/AdA 



fx AdA 

versus the wavelength-based asymmetry statistic: 

(Aio,r Ap) 

(Ap — A 10j b) 



(3) 



(4) 



for our sample, where A p is the wavelength of the peak of the emission line, and A 10) & and A 10 r 
are the wavelengths at which the line flux first exce eds 10% of the peak on th e blue side and 
on the red side of the emission line, respectively (see Rhoads et al. 2003 . 2004 . and Paper 
Each of the confirmed Lya emitters in this sample satisfies af > 1.0 or a\ > 1.0, and 52 out 
of 59 sources satisfy both. As we found for the lower resolution Keck/LRIS sample in Paper 
I, the present Lya sample, observed with higher spectral resolution using Keck/DEIMOS, is 
systematically segregated from low-redshift [O II] A3727 in df-a\ space. 

As a final diagnostic, we note that in each of our confirmed Lya-emitters for which the 
continuum is sufficiently well-detected, the spectrum shows a continuum decrement consis- 
tent with the onset of absorption by the Lya forest at A rcst = 1216 A. The break amplitude 
is typically characterized by 1 — /* ort / 'i^ 0118 , where we define /* hoIt as the variance-weighted 
flux density in a 1200 A window beginning 30 A below the emission line; f l ° ns is the same, 
but above the emission line. In the 24 sources for which /^ ong is detected to better than 2a, 
all but two sources have 1 — /* ort / P° ns > 0.5, consistent with continuum break amplitudes at 



5 As in Paper I, the error bars on a\ and a/ were determined with Monte Carlo simulations in which 
we modeled e ach e mission line with the truncated Gaussian profile described in iHu et al.l (|2004) and 



Rhoads et al.l (|2004f ). added random noise in each pixel according to the photon counting errors, and mea- 



sured the widths a(a\) and cr(a/) of the resulting distributions of a\ and a/ for the given line. That is, for 
each a\ t i, the error 5a\^ = u(aA.i), and similarly for each a/^. 
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4.5 in theoretical models (e.g.. lMadaulll995l ; IZhang et al.lll997l ). in the lower resolution 



Keck /LRIS sample presented in Paper I, and in other similar datasets (see lStern &: Spinrad 



19991 . and references therein). 



4. Discussion 



Together with the observations presented in Paper I (and accounting for minor overlap 
in the samples), the 59 Lya-emitters confirmed herein bring the total catalog of spectroscop- 
ically confirmed z ~ 4.5 Lya-emitting galaxies to 73 objects in the ~ 0.7 degrees 2 surveyed 
by the LALA imaging. We now update the characteristics of this population as they were 
estimated in Paper I by investigating the distribution of the total sample in equivalent width. 
We then construct a z ~ 4.5 Lya luminosity function, carefully accounting for survey in- 
completeness and for spectroscopic sensitivity, and we compare the result to Lya luminosity 
functions spanning 3.1 < z < 6.6. 



4.1. The Equivalent Width Distribution 

As in Paper I, we determine the rest-frame equivalent widths directly from the spectra 
according to W^ st = {Fgj fx,r) / '(1 + z), where Ft is the flux in the emission line and f\^ r is 
the measured red-side continuum flux density. The resulting equivalent width distribution 
is plotted in Figure [3, together with the equivalent widths measured in Paper I. 

Before interpreting this distribution, one should be cautioned that the W{ cst determi- 
nation is very sensitive to uncertainty in the measured continuum. Since the continuum 
estimate enters into the denominator of the expression for W{ cst , the characteristically small 
continuum values and their large fractional uncertainties cause significant scatter in the 
measurement, and the resulting error is neither Gaussian nor symmetric about the mea- 
sured value. Especially problematic is the fact that the largest values of W{ cst are also the 
least certain. Detailed discussions of the uncertainties in measuring W{ cst in high-redshift 
Lya-emitters, along with the co mplicating effect s of dust content, gas kinematics, and inter- 



galactic absorption, are given in iHu et al.1 (120041 ) and in Paper I. 



With these caveats in mind, we rigorously treated the error bars on the equivalent width 
estimates, and we restricted the analysis to sources with red-side continuum signal-to-noise 
ratios > 1. To determine the equivalent width error bars, we first associated each measured 
line flux ± 5Fi ti with a Gaussian probability density function (PDF) centered on F^ 
with width o = 8Fn; we proceeded similarly for the measured continuum fluxes. We then 
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generated a grid of line flux versus continuum flux on which each node has an associated 
equivalent width and is assigned a weight given by the probability distribution on each of 
its flux axes. Next we collapsed the grid into a histogram of equivalent widths, adding 
the weight from each grid point to the appropriate equivalent width bin. The result is a 
non-Gaussian PDF Pi(w) for which Pi(w) dw is the probability of observing W{ c f in the 
interval w < W™f < w + dw. The error bars 5w + and 5w_ are then la confidence intervals 
determined by integrating over the probability density functions Pi(w). They are symmetric 
in probability density-space in the sense that f w _ g Pi{w') dw' = f™ +Sw+ Piiw') dw' = 0.34. 

We find the resulting distr ibution to be broad ly consistent wit h the equivalent widths 
presented in iFujita et al.l (120031 ) for z ~ 3.7 and in iHu et al.l (120041 ) for z ~ 5.7. While the 
majority of sources can be understood as comparatively young (1 to 10 Myr) galaxies with 
Salpeter initial mass functions (IMFs), a non-negligible fracti on exceeds the largest rest- 
frame equivalent widths expected from such stellar populations. iMalhotra &: Rhoadd (120021 ) 
use a Salpeter initial mass function, an upper mass cutoff of 120 M , and a metallicity 
of l/20th solar to find maximum Lya equivalent widths of 300 A, 150 A, and 100 A for 
stellar populations of ages 10 6 , 10 7 , and 10 8 years, respectively. Adopting a correction factor 
of 0.64 as an upper limit to the effect of IGM absorption on the measurement of W{ cst in 
spectroscopy effectively reduces these upper limits to 190 A, 100 A, and 60 A( see discussion 
in Paper I). Owing to the lower me tallicity used in their models, the pre-IGM-corrected 
values of IMalhotra fc Rhoadd (120021 ) are slig htly higher than the canonical limiting Lya 
rest-frame equivalent width of 240 A given by I Char lot fc Falll ( 1993 ). 



Using the ensemble of Pi{w) described above, we find that 12% - 27% (90% confidence) 
of the galaxies in this sample show W{ cst > 240 A, and 17% - 31% (93% confidence) show 
W^ est > 190 A. Both results are nearly identical to the values given in Paper I. On the 
simplest interpretation, these galaxies are required to be very young (age < 10 6 years), or 
to have IMFs skewed in favor of the production of massive stars. The possibility that AGNs 
in our sample are producing stronger-than-expected Lya emission seems unlikely due to 
the comparatively narrow velocity widths of the Lya lines and to the absence of the high- 
ionization state UV emission lines symptomatic of AGN activity. Mor eover, deep (^170 ks) 
Chandra/ ACI S imaging of L AL A z « 4.5 candidates in both Bootes (IMalhotra et al.l 120031 ) 



and in Cetus (jWang et al.ll2004l ) resulted in X-ray non-detections to an average 3a limiting 



luminosity of L2_8kev < 2.8 x 10 erg s . This limit is roughly an order of magnitude 



Stern et al. 


2002; 


Norman et al. 


2002; 



on 



the typical X-ray to Lya luminosity ratio for the L ya galaxy sampl e to th e obs erved values 



of this ratio for quasar and Seyfert galaxy samples, IMalhotra et al.l (120031 ) and IWang et al. 
(120041 ) conclude that AGN account for < 5% of the Lya galaxy sample. 
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4.2. Empirical Cumulative Luminosity Function 



In Figure we present an empirical cumulative Lya line luminosity function computed 
for our sample at z ~ 4.5 and compare this to luminosity functions computed for several 
other samples spanning 3.1 < z < 6.6. The cumulative luminosity function gives for each 
Lya line luminosity L(Lya) the total number density of Lya lines brighter than L(Lya). 
The comparison samples are drawn from spectroscopic follow-up of narrow-band surveys with 
roughly comparable flux limits and candidate-selection criteria (except where n oted, below). 



We d o not include non-spectroscopic Lya-emitter luminosity functions (e.g., lOuchi et al. 



20031 ) among the comparison samples. In each case, we converted the reported Lya line 
fluxes to line luminosities using a A-cosmology with Q M = 0.3 and = 0.7, and H = 70 
km s _1 Mpc -1 , and we made a minimal attempt to account for incompleteness^. Specifically, 
the volume from which Lya-emitting candidates were selected by their narrow-band excess 
is simply defined by the solid angle covered by the narrow-band imaging and the redshift 
range allowed by the narrow-band filter. However, the effective volume surveyed by the 
spectroscopic follow-up is smaller than the imaging survey volume by a factor of N spcc / N ca _ nd , 
where iV cand is the total number of Lya-emitting candidates discovered in the imaging, and 
iV" spec is the number of candidates actually targeted for spectroscopy. We estimated the un- 
certainties in the cumulative luminosity functions with Monte Carlo simulations. Assuming 
the errors in the Lya line fluxes are Gaussian, we created synthetic data sets by drawing 
randomly from the Gaussian Lya flux PDFs for each object in each sample. The result for 
each sample was then a distribution of cumulative luminosity functions, which may be used 
to define upper and lower confidence intervals. Figure [8] depicts 95% confidence intervals; 
where more than one survey is plotted, just the confidence intervals for the survey with the 
largest range in line fluxes is depicted. 

No strong evolution is readily evident in the cumulative Lya luminosity functions be- 
tween z ~ 3 and z ~ 6. The only significant scatter between luminosity functions occurs 
between the various z ~ 3 surveys, and that scatter likely finds its origin in differences 
in t he manner in which the experiments were performed. Foremost, the area surveyed by 
the ICowie fc Hul (119981 ) effort is comparatively sm all: just 25 arcm i n 2 in each of two fields 
( HDF and SSA 22), as opposed to 300 ar cmin 2 in iKudritzki et al.l (120001 ) and 132 arcmin 2 



m 



Fujita et al.l (120031 ) . ICowie fc Hul (119981 ) note that the number counts in their HDF field 



appears to be 2.5 times richer in narrow-band excess objects than their SSA 22 field, high- 
lighting the susceptibility of small survey areas to cosmic variance. Separately, as noted by 



Hu et al.1 (|2004l ) provide fluxes for their 



5.7 sources as measured in narrow-band ima ging;, rather 
than L yq line fluxes as measured in spectroscopy. As such, we adopt the conversion given in Istern et al.l 
|2005l ) to estimate / LyQ from /^ B . 
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Hu et al.l (120041 ). the iFujita et al.l (120031 ) data may comparatively under-represent the den- 



sity of Lya-emitters due to their m ore stringent equiva lent width criterion of W^ hs > 250 



A, as opposed t o W? bs > 77 A in ICowie fc Hul (119981 ) and effectively W{ hs > 100 A in 



Kudritzki et al.l (j2000h . 



4.3. The V/V^ Estimate 

We now perform a more rigorous measurement of the z ~ 4.5 Lyq luminosity fun ction 
using a modified version of the V/V max method (e.g., Hogg et al. 1998 : Fan et allboOll ). For 



each galaxy, V msac is the volume over which Lyct of a given luminosity could be located and 
still be detected by our survey; the luminosity function is then the sum of the inverse volumes 
of all galaxies in the given luminosity bins. Our modifications to the V/V m ^ method account 
for incompleteness in two senses. First, not e very galaxy c andid ate identified in imaging was 



targeted in follow-up spectroscopy. Following iHogg et al.l (119981 ). Figure [HI shows the fraction 



of narrow-band selected candidate Lya-emitters which were targeted for spectroscopy as a 
function of flux in the band in which the candidate was detected. We label this a priori 
completeness function r/ try ; the candidate Lya flux f hya can be roughly estimated from the 
flux in the narrow-band f™ B with f hya = w n (f™ B — f*), where w n is the width of the narrow- 
band filter and is the flux of the candidate in the .R-band. 

Second, even if a candidate Lya-emitter was selected for spectroscopy, its inclusion in 
the luminosity function depends on the detection and identification of the Lya line. Our 
spectroscopic sensitivity to Lya emission as a function of flux and redshift is shown in 
Figure HI we label this function p dctcct . As discussed in section § 13. 1\ p dctcct can be interpreted 
as the probability that a putative Lya emission line of a given flux and a given redshift 
would have been detected in our spectroscopic campaign. 

In the presence of these selection effects, the available volume for a galaxy with Lya 
emission of flux f Lya is 

f Z2 d 2 V c 

Knax = J Vtry(fl ya ) P d , tc Mly a , A ^^7^ ^ ' ( 5 ) 

where the comoving volume element in a solid angle dfl and redshift interval dz is the familiar 

d 2 v r ( c \ 3 r r dz' i 2 i 

(6) 



dVtdz \H n ) |./n E(z') I E(z) 

with 



e(z) = {n u (i + zf + n k (i + zf + n A } 1/2 . (7) 
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In equation [51 AQ is the solid angle covered by the LALA survey, and /' is the Lya line 
flux for the source in question if it were located at redshift z'. The lower limit of integration 
Z\ is set by the lowest wavelength at which Lya could be detected by our narrow-band filters, 
corresponding to z ~ 4.37. The upper limit of integration z 2 is set in one of two ways. If the 
Lya luminosity for a source is bright enough that the line remains above the survey flux limit 
out to the highest redshift accessible by our filter set, then z 2 is simply equal to the upper 
redshift limit for the survey, z ~ 4.57. For fainter sources, z 2 is taken to be the redshift at 
which the Lya flux falls below the survey flux limit; in this case, 4.37 < z 2 < 4.57. 

Having computed V maK for each galaxy, we may compute the differential Lya luminosity 
function $(L), the number density of galaxies per logarithmic interval in Lya luminosity. In 
a given luminosity bin of width AlogL centered on Lj, this is given by 



AlogL 



E 



Here, the index i denotes the luminosity bin and j denotes the galaxies within the bin, where 
the galaxies summed in a given bin are selected by their Lya luminosities according to 



log Lj - log ^ | < 



AlogL 



(9) 



Finally, the uncertainty in the luminosity function may be estimated with 

1/2 



AloeL 



E 



(10) 



In Figure QUI we present the resulting differential Lya luminosity function at z ~ 4.5. 
We also fit the data with a Schechter function. If $(L) dL is the comoving number density 
of galaxies with luminosities in the range (L, L + dL), then the corresponding Schechter 
function is 

W^=^(^) Q exp(-^) dL, (11) 

where $* is the normalization, L* is the characteristic break luminosity, and a sets the slope 
at the faint end. This is related to the number density of galaxies in logarithmic intervals by 

* (i)<i(iosi) = Gdb) (£) (f)° cxp (4) d(logL) • 



and it is this function which we fit to our data. As in Ivan Breukelen et al.l (120051 ). because 
the binned data points are few, we choose to fix a = —1.6 so as to fit with only two free 
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parameters, $* and L*. Thi s choice fits well with the luminosity distribution of both LBGs 
and Lya-emitters at z ~ 3 (jSteidel et al.lll999l . l2000l ). We find best-fit luminosity function 
paramters L* = (10.9 ±3.3) x 10 42 erg s" 1 and $* = (1.7±0.2) x 10" 4 Mpc~3 (or equivalently, 
log(L*) = 43.04 ± 0.14 and log($*) = -3.77 ± 0.05). The error bars on L* and $* are the 
la formal errors computed from the covariance matrix in the nonlinear least-squares fit, 
scaled by the mea sured value of x 2 ■ That is, 5L* = ol* a/ X 2 / n DOF, and similarly for 5$* 
rtPress et al.lll992h . 



Our z ~ 4.5 sample provides one of the best measured Lya luminosity functions to date. 
We can study redshift evolution of the Lya luminosity function by comparing to results from 
the literature. Recognizing that the uncertainties in L* and $* are strongly correlated, we 
examine not only the individual parameters but also the product which is proportional 

to Lya luminosity density, and which generally has smaller uncertainties than the individual 
parameters. For our sample, log(L*$*) = 39.27. 



At lower redshift, there is a z ~ 3.1 LF by (IGronwall et al.l 120071 ). who fit all three 



parameters. They find a = -lA9t° H, log(L*) = 42M±° f 5 , and log($*) « -2.89 ± 0.04, 



whence log(L*$*) ^ 39.7 5. A t the high redshift end , we compare to LFs at z = 6.5 by 



Malhotra fc Rhoads (2004) and 



Kashikawa et al.l (120061 ). and at z = 5.7 by lMalhotra fc Rhoads 



( 120041 ) and IShimasaku et al.l (120061 ) , all derived by fixing the faint end slope a = —1.5 and 



fitting L* and $*. At z — 6.5, the LFs are similar to our z 
(|2004h find log(L*) = 42.6, log($*) = -3.3, and log(L*$ 



4.5 result: Ma 



hotra fe Rhoads 



Kashikawa et al. 



= 39.3, while 

( 120061 ) find (for their combined spectroscopic + photometric samp l e) lo g(L*) = 42.6, log($*) = 
-2.88, and log(L*$*) = 39.72. At z = 5.7. iMalhotra fc Rhoadsl toOf\ fin d log(L*) = 43.0, 
log($*) = -4.0, and log(L*$*) = 39.0, while IShimasaku et al.l J2006h find log(L*) = 
42.9 ± 0.14, log($*) = -3.2 ± 0.17, and log(L*$*) = 39.7. The obvious differences be- 
tween the LFs at each redshift may be caused by any combination of (a) simple uncertainty 
in deriving the LF from modest sized samples; (b) field-to-field variations in Lya galaxy 
density; or (c) differences in the methods used to derive Sc hechter function p aram eters, 
and in part to l o cal y ariations in Lya galaxy density. The Kashikawa et al.l (120061 ) and 
(2006 1) LFs are der i ved fr om larger total samples, but from a single survey 



Shimasaku et al. 



field, while the IMalhotra fc Rhoadsl (120041 ) LFs are based on a combination of several older, 
smaller samples from a few widely separated fields. Regardless, if we take the difference 
between these various z ~ 6 LFs as an empirical indication of total present uncertainties, 
the z ~ 4.5 LF derived in the present paper supports a roughly constant Lya luminosity 
density over the range z = 4.5 ± 1.5. 
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4.4. Comparison to LBGs 



It is interesting to compare this result to the evolution of the rest-UV luminosity density 
and cosmic star formation rate density (SFRD) derived from LBGs over the same redshift 
range. Estimates of the z ~ 6 SFRD based on the Great Observatories O rigins Deep Sur- 
vey/Advanced Ca mera for Surveys (G OODS/ACS; iGiavalisco et al.ll2004al ) and the Hubble 
Ultra Deep Field (IBeckwith et al.ll2006l) show a factor of 1-5-6 dro p between z ~ 3 and z ~ 6 
(feunker et allbooi feouwens et al.lbood : IGiavalisco et"al"1l2004bh . 



What does this mean for the Lya-selected galaxies? ISteidel et al.l (120001 ) report that 
Lya-selection with an equivalent width criterion typical of narrow-band surveys would return 
20% - 25% of their z ~ 3 LBGs. If Lya-emitters are merely a subset of the LBG population 
which happen to have been detected during a stage of strong Lya production, than we 
would expect the Lya luminosity density to decline beyond z ~ 3, in step with the global 
star formation rate density. Integrating the luminosity functions discussed in section 14.31 
shows no compelling evidence for such a decline. Though a modest decline cannot be firmly 
ruled out, we may nontheless speculate that the Lya-emitters as a population are evolving 
differently from the LBGs. 

The bolometric luminosities of Lya galaxies are typically lower than Lyman break galax- 
i es, and provide a hint that they are less mass i ve. Detailed spectr al energy distribution fitting 
( IGawiser et all 120061 ; iFinkelstein et al.l 120071 ; iPirzkal et al.l 120071 ) bears out this preliminary 
inference, showing typical masses of ~ 1O 8 M and a ges ~ 10 7 to 10 8 years. The correla - 
tion strengths of Lya galaxies and LBGs are similar (lOuchi et al.ll2003l ; IKovac et al.l 120071 ). 
indicating similar masses of haloes. From the expected halo mass one can predict volume 
number density of Lya emitters. Comparing the expected an d observed numbe r densities 
implies a duty cycle of Lya-emission in the range 6% - 50% (IKovac et al.l 120071 ) . A simi- 
lar duty cycle, 7.5% - 15%, is infe rred from stellar population modelling of the photomeric 



sample flMalhotra fc Rhoadsll2002l ). 



4.5. Implications for Reionization 



The spectroscopic observations of the z > 6 quasars yielded the first detections of the 
l ong-awaited Gunn - Peterson trough, implying at least the end of reionization at z rs 6 
( iBecker et al.l 1200 ll ; iDjorgovski et al.l l200ll ; iFan et al.ll2002l ). Subsequently, the Wilkinson 
Microwave Anisotropy Probe (WMAP) identified a large amplitude signal in the temperature- 
polarization maps of the cosmic microwave background (jSpergel et al.ll2003l ; iPage et al.ll2006l ) 
indicating a large optical depth to Thomson scattering and favoring reionization instead at 
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z ~ 11. The WMAP results are not necessarily inconsistent with those of the quasar Gunn- 
Peterson troughs. Only a small neutral fraction (x^ M ~ 0.001) is required to produce 
the Gunn- Peterson effect, so one plausible scenario is that reionization may have been an 
extended event, beginning early but not completing until z ~ 6. Alternatively, a variety of 
theoretical models now suggest that reionization occurred twice, first at z ~ 20 with the onset 
of zero-metallicity Population III s tars, and then again by massive Population II stars forme d 
after a partial recombination (e.g. , ICenll2003l ; lHaiman fc Holderll2003l ; ISomerville et al.ll2003l ) . 



High-redshift Lya-emitting galaxies offer another perspective on this issue, as the vis- 



Haiman & Spaans 


1999; 


Santos 


2004 


)■ 



present first attempts to exploit this fact by comparing luminosity functions of Lya-emitters 
at z ~ 5.7 and z ~ 6.6. They find no measurable evolution between these epochs, from which 
the y infer that the IGM re mains largely reionized from the local universe out to z m 6.5 (but 



sec 



Haiman fc Cenll2005l ). iKashikawa et al.l (120061 ) . applying the same test, find possible evi- 
dence for observed Lya LF differences between z = 5.7 and 6.5 at the fact o r of 2 level. They 
suggest neutral gas at z ~ 6.5 as the explanation, though iDijkstra et all (120071 ) argue that 



the observations could equally well be explained by the ongoing growth of cosmic structure 
from z = 6.5 to z = 5.7. 



By using the Lya galaxy sample from iTaniguchi et al.l (120051 ). iMalhotra fc Rhoads 

( 120061 ) showed that at least 30% of the IGM by v olume is ionized at z « 6.5. This is 
corroborated by dark gap statistics in GP troughs jFan et alJlioO^ I. All the Lya tests of 
reionization assume that there is no intrinsic evolution in the Lya luminosity functions be- 
tween z = 5.7 and 6.5. In this paper we show that there is little evolution in Lya luminosity 
function from z = 6.6 to z = 3.1, thus strengthening the conclusion that the IGM is not 
substantially neutral at z — 6.5. 

Significantly, the related question of what is responsible for reionization remains at 
large. It has long been recognized that AGN at earl y epochs are insufficient, owing to thei r 
rapid decline in space density at high redshift fe.g.. iMadau et al.lll999l ; iBarger et al.ll2003l ). 
Based on their analysis of the UDF, iBunker et al.l (12004 ) conclude that the cosmic SFR 
in directly observed z ~ 6 L BGs was roughly five times too low to reionize the Universe. 
Yan fc Windhorstl (120041 ) and iBouwens et all (120061 ) argue that the ionizing photon budget 
is sufficient provided one accounts for sample incompleteness using a sufficiently steep slope 
at the faint end of the luminosity function. IMalhotra et al.l (120051 ) argue that the ionizing 
flux density may be very inhomogeneous due to large scale structure, as seen in galaxies in 
the Hubble Ultra Deep Field, and that the directly observed galaxies at z ~ 6 do produce 
sufficient photons for reionization in overdense regions. 



-18- 



We estimate that the contribution to the cosmic SFR from Lya-emitters at this epoch 
is lower than that of the LBGs (p SFR (Lya) ~ 0.003 M yr^ 1 Mpc -3 , as compared to 
p SFR (LBG) ~ 0.005 M Q yr^ 1 Mpc~ 3 ) when integrated over the same limits. Consequently, 
though high-redshift Lya-emitters are proving to be a useful probe of the history of reion- 
izatio n, they are evident ly not its cause. While extinction corrections could occasionally be 
large (jChary et al.ll2005l ) and could modify this conclusion, most well st udied Lya galaxies 
have very modest extinction (IFinkelstein et al.ll2007l ; iPirzkal et al.ll2007l ). 



When we compare the luminosity function of Lya-emitters at z ~ 4.5 to luminosity 
functions for similarly assembled samples spanning 3.1 < z < 6.6, we find no evidence 
for ev oluti on over thes e epoc hs. This result bolsters the conclusion by iMalhotra fc Rhoads 
(120041 ) and lStern et al.l (120051 ) that the IGM remains largely reionized from the local universe 
out to z ~ 6.5. However, it is somewhat at o dds with the factor of 1.5-6 drop in the 



cosmic star formation rate density measured by (IBunker et al.l 120041 ; iBouwens et al.l 12006 
Giavalisco et al.l l2004al ) between z ~ 3 and z ~ 6 in Lyman- break galaxies selected in the 
exceptional imaging of the UDF. It seems that these two populations — Lya emitters and 
Lyman Break Galaxies — follow different evolutionary histories. The disentanglement of this 
issue will likely rely on extensive followup observations of large samples, so that we can study 
the continuum and absorption lines of many Lya galaxies, and conversely the Lya properties 
of the break-selected galaxies. 
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Fig. 1. — Distribution of redshifts for spectroscopically confirmed Lya emission lines in the 
Cetus field (left; 02:05:20 -04:55, J2000.0) and in the Bootes field (right; 14:25:57 +35:32, 
J2000.0). The redshifts labeled "DEIMOS" denote galaxies confirmed with our campaign of 
Keck/DEIMOS spectroscopy, described in this paper. The redshifts labeled "LRIS" denote 
galaxies confirmed with our campaign of Keck/LRIS spectroscopy, described in Paper I. The 
overlays are arbitrarily scaled transmission curves for the five narrow-band filters employed 
in the imaging component of this survey. 
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Fig. 2. — Sample spectra from the set of 59 z xs 4.5 Lya-emitting galaxies confirmed with Keck/DEIMOS, 
with a wavelength range selected to highlight the emission-line profile. The measured redshifts and asym- 
metry statistics (§ I3.2j) are indicated in the upper right of each panel. The representative error bar (upper 
left) is the median of the flux error in each pixel over the wavelength range displayed. The spectra have been 
smoothed with a 3-pixel boxcar average. 
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Fig. 3. — Spectroscopic success rate as a function of the flux in the narrow-band in which the 
candidate was selected. This plot combines the results of the Keck/DEIMOS observations 
made for this paper and the Keck/LRIS observations described in Paper I. 
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Fig. 4. — Empirical, cumulative distribution of spectroscopic sensitivity to Lja emission, as 
a function of source redshift and Lya flux. The contours span 10% to 90% in 10% steps. The 
dark lines denote the 10%, 50%, and 90% contours. The distribution is plotted cumulatively 
so that it can be interpreted as the probability that a putative Lya emission line of a given 
flux and redshift would have been detected in our Keck/DEIMOS spectroscopic campaign. 
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Fig. 5. — Sample Lya emission line profile (left) compared to two common low-redshift in- 
terlopers: [O II] A3727 (center) and [O III] A5007 (right). The top figure in each case is the 
one-dimensional extracted spectrum; the bottom figure is a section of the two-dimensional 
data from which it was extracted. Note that we resolve the [O II] A3727 doublet with 
our Keck/DEIMOS spectroscopic setup, thereby eliminating [O II] A3727 as the main low- 
redshift interloper in our survey. The [O III] A5007 line can typically be identified by neigh- 
boring [O III] A4959 at one-third its strength, or by neighboring B./3. 
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Fig. 6. — Scatter plot comparing the flux-based asymmetry statistic a/ and the wavelength-based asymme- 
try statistic a\ of known high-redshift Lya-emitters to a sample of [O II] A3727-emitters at z ~ 1, updated 
from Paper I. The points labeled "DEIMOS" denote galaxies confirmed with our campaign of Keck/DEIMOS 
spectroscopy, described in this paper. The points labeled "LRIS" denote galaxies confirmed with our cam- 
paign of 400l/mm-grating Kec k/LRIS spectroscopy, described in Paper I. The thr ee Lya-emitter s at z = 5.7 



are from 



Stern et al 



Rho ads et al.l (|2003l ), and the two Lya-emitters at z = 6.5 are from iRhoads et al.l ( 2004) and 



(2005). The 28 [O II] A3727-emitters at z ~ 1 were provided by the DEEP2 team (jPavis et al 



200.4 A. Coil 2004, private communication); their Keck/DEIMOS 1200^/mm- grating spectra were smoothed 



to the Keck/LRIS 400£/mm-grating resolution by convolution with a Gaussian kernel. The representative 
error bar (lower right) is the median of the errors on the individual a/ and a\ for the combined Keck/LRIS 
and Keck/DEIMOS sample. 
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Fig. 7. — Histogram of the spectroscopic rest-frame equivalent widths for the z = 4.5 pop- 
ulation, determined with W^ cst = (Fg/ ' f\ <r ) / '(1 + z), where Fg is the flux in the emission line 
and f\ tr is the measured red-side continuum flux density. The sources labeled "DEIMOS" 
denote galaxies confirmed with our campaign of Keck/DEIMOS spectroscopy, described in 
this paper. The sources labeled "LRIS" denote galaxies confirmed with our campaign of 
Keck/LRIS spectroscopy, described in Paper I. Representative error bars on the equivalent 
widths are plotted at left and at right. Notably, the highest equivalent widths are generally 
the least certain, as they correspond to the faintest (and hence least certain) continuum 
estimates. 
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Fig. 8. — Comparison of empirical, cumulative Lya luminosity functions computed with only minimal 
completeness correction for several spectroscopic surveys spanning 3.1 < z < 6.6. The cumulative luminosity 
function gives for each Lya line luminosity L(Lya) the total number density of Lya lines brighter than 
L(Lya). The shaded regions represent 95% confidence intervals based on the Monte Carlo simulations 
described in § 14.21 Where more than one survey is plotted, just the confidence intervals for the survey with 
the largest range in line fluxes is depicted. No strong evolution is evident over the redshift range depicted. 
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Fig. 9. — Probability as a function of narrow-band flux that a candidate Lya-emitter was 
targeted for spectroscopy, divided into the fraction of targets that were confirmed and the 
fraction of targets that were not. 
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Fig. 10. — Differential Lya luminosity function for our z = 4.5 sample, computed using 
the V/V maiK method. The sample includes both the Keck/DEIMOS data presented in this 
paper and the Keck/LRIS data presented in Paper I. The error bars are the la statistical 
uncertainties given by the root variance shown in equation [10J The background histogram 
(dashed) gives the number of individual sources contributing to each luminosity bin. The 
solid curve shows the best fitting Schechter function model, with L* = (10.9 ± 3.3) x 10 42 
erg s -1 , $* = (1.7 ± 0.2) x 10" 4 Mpc~3, and a fixed faint-end slope a = —1.6. 
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Table I. Spectroscopic Properties 
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4-30 


7.4 


± 


0.6 


260 


0.011 ± 0.046 


0.126 


± 


0.053 


J142525 


.9+352349 


4 


.471 


3.27 


± 


1.34 


33±} 6 3 


7.0 


± 


0.7 


240 


0.072 ± 0.048 


0.267 


± 


0.057 


J142526 


.2+352531 


4 


.464 


2.76 


± 


1.13 


r, -4-105 
85 l39 


6.2 


± 


0.4 


190 


0.067 ± 0.050 


0.087 


± 


0.054 


J142531 


.8+352652 


4 


.482 


0.94 


± 


0.40 


-1 0+9 


7.4 


± 


1.5 


270 


0.034 ± 0.035 


0.140 


± 


0.041 


J142532 


.9+353013 


4 


.534 


5.49 


± 


1.00 


20ll™ 


7.1 


± 


0.3 


250 


0.005 ± 0.015 


0.075 


± 


0.018 


J142535 


.2+352743 


4 


.449 


6.23 


± 


2.54 


159±™ 


6.0 


± 


0.2 


180 


0.001 ± 0.048 


0.106 


± 


0.057 


J142539 


.5+353902 


4 


.432 


1.52 


± 


0.67 


240 +2182 
zw -126 


4.0 


± 


1.6 


< 200 h 


0.049 ± 0.019 


0.017 


± 


0.022 


J142541 


.7+353351 


4 


.409 


3.24 


± 


1.34 


108l« 9 


5.4 


± 


0.8 


130 


-0.042 ± 0.045 


0.080 


± 


0.050 


J142542 


.0+352557 


4 


.393 


1.05 


± 


0.44 


30±T 3 


6.4 


± 


1.4 


210 


0.028 ± 0.033 


0.092 


± 


0.042 


J142542 


.6+352626 


4 


.450 


1.49 


± 


0.62 


19 +1 7 ° 


7.5 


± 


0.9 


270 


0.101 ± 0.043 


0.215 


± 


0.059 


J142544 


.5+354325 


4 


.533 


2.84 


± 


1.20 


1,1+129 
1J1 -60 


7.4 


± 


1.1 


260 


0.002 ± 0.018 


0.059 


± 


0.030 


J142546 


.8+354315 


4 


.443 


0.72 


± 


0.33 


^ u -19 


6.3 


± 


1.5 


200 


0.012 ± 0.016 


0.049 


± 


0.022 


J142547 


.8+354200 


4 


.539 


1.11 


± 


0.48 


> 56 g 


4.9 


± 


0.8 


90 


0.029 ± 0.017 


0.007 


± 


0.024 


J142548 


.4+352740 


4 


.546 


1.21 


± 


0.50 


> 24 g 


5.4 


± 


0.7 


130 


-0.058 ± 0.038 


0.004 


± 


0.067 


J142555 


.4+353039 


4 


.423 


10.31 ± 1.88 


560+ 467 


10.8 ± 0.6 


450 


0.025 ± 0.015 


0.049 


± 


0.023 


J142556 


.7+354234 


4 


.425 


2.26 ± 0.96 


189+ 284 
162+lf 


4.8 ± 0.6 


80 


-0.019 ± 0.020 


0.032 


± 


0.022 


J142556 


.8+354215 


4 


,426 


2.85 


± 1.20 


5.8 


± 0.7 


170 


0.021 ± 0.020 


0.047 


± 


0.022 


J142559 


.8+353513 


4 


.394 


1.39 ± 0.28 


158ljg 9 


5.6 ± 0.6 


160 


0.005 ± 0.013 


0.023 


± 


0.014 
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Lyo Flux b 


■pyrestc 


FWHM d 


Av c 


Continuum (/ijy) f 


Continuum (/^Jy) f 


Target 






(10~ 17 erg cm~ 2 s _1 ) 


(A) 


(A) 


(km s-i) 


Blue Side 


Red Side 


11 42559 8+353748 


4 


.420 


4.95 


± 


0.91 


55+10 


7.0 ± 0.5 


240 


109 + 018 


240 + 0019 


J142601.3+353618 


4 


.475 


1.41 


± 


0.27 


> 171 s 


5.0 ± 0.4 


100 


-0.013 ± 0.011 


-0.003 ± 0.013 


J142602.0+354554 


4 


.473 


1.83 


± 


0.78 


85 -38 


4.8 ± 0.7 


70 


-0.003 ± 0.019 


0.058 ± 0.025 


J142612. 2+353541 


4 


.418 


1.90 


± 


0.36 


> 140 s 


6.3 ± 0.5 


200 


0.082 ± 0.013 


0.002 ± 0.017 


J142624.4+353832 


4 


.460 


2.46 


± 


0.46 


OQfi+1182 
JZU _145 


5.4 ± 0.5 


140 


-0.012 ± 0.017 


0.021 ± 0.022 


J142627.5+353717 


4 


.488 


2.11 


± 


0.43 


4 o+24 
4Z -14 


6.5 ± 2.4 


210 


0.044 ± 0.025 


0.135 ± 0.051 


J142628.5+353809 


4 


.409 


3.82 


± 


0.71 


> 143 s 


6.4 ± 0.4 


210 


-0.038 ± 0.029 


-0.011 ± 0.041 


J142653.5+353356 


4 


.494 


2.12 


± 


0.77 




7.8 ± 1.8 


280 


-0.016 ± 0.015 


0.040 ± 0.017 


J142658.8+353144 


4 


.495 


2.16 


± 


0.79 


31+11 


6.3 ± 1.0 


200 


-0.019 ± 0.015 


0.191 ± 0.018 


J142706.3+353224 


4 


.480 


0.99 


± 


0.38 


ioo+if 


5.0 ± 1.6 


100 


0.021 ± 0.014 


0.027 ± 0.018 


J142709.1+352738 


4 


.407 


1.77 


± 


0.65 




8.4 ± 1.4 


320 


0.039 ± 0.014 


0.056 ± 0.017 


J142709.2+352409 


4 


.520 


1.62 


± 


0.59 


> 85 s 


5.8 ± 1.2 


170 


0.029 ± 0.015 


0.008 ± 0.022 


J142709.8+352641 


4 


.405 


1.78 


± 


0.66 


82±|| 


6.0 ± 0.8 


180 


0.023 ± 0.017 


0.027 ± 0.021 


J142712. 2+353029 


4 


.380 


2.35 


± 


0.86 


4.3 ± 2.4 


< 200 h 


0.035 ± 0.017 


0.076 ± 0.021 



a The redshift was derived from the wavelength of the peak pixel in the line profile smoothed with a 3-pixel boxcar average. We 
estimate the error in this measurement to be <5 Z ^ 0.0005, based on Monte Carlo simulations in which we added random noise to each 
pixel of every spectrum according to the photon counting statistics, and then re-measured the redshift in each case. This measurement 
may overestimate the true redshift of the system since the blue wing of the Lyo emission is absorbed by foreground neutral hydrogen. 

b The line flux was determined by totaling the flux of the pixels that fall within the line profile. No attempt was made to model 
the emission line or to account for the very minor contribution of the continuum to the line. Quoted uncertainties account for photon 
counting errors alone, excluding possible systematic errors. Despite these caveats, the Lya line fluxes measured from the spectra agree 
with narrow band imaging to la in all but three cases. 

c The rest frame equivalent widths were determined with W^ ast = //>, r )/(l + z), where Ft is the flux in the emission line and f\ r 
is the measured red-side continuum flux density. The error bars 8w+ and 6w- are la confidence intervals determined by integrating 
over the probability density functions Pi(w) described in £j 14.11 The error bars are symmetric in probability density-space in the sense 
that £_ Sw _ P(w>) dw' = f™ +Sw + Pi(w') dw> = 0.34. 

d The FWHM was measured directly from the emission line by counting the number of pixels in the unsmoothed spectrum which 
exceed a flux equal to half the flux in the peak pixel. No attempt was made to account for the minor contribution of the continuum 
to the height of the peak pixel. The error bars were determined with Monte Carlo s imula tions in which we modeled each emission line 
with the truncated Gaussian profile described in iHu et alj d2004ll and lRhoads et al.l l l2004h . added random noise in each pixel according 
to the photon counting errors, and then measured the widths cr(FWHM) of the resulting distribution of FWHM for the given line. 

°The velocity width Av was determined by subtracting in quadrature the effective instrumental resolution for a point source, and is 
therefore an upper limit, as the target may have angular size comparable to the i; 1" seeing of these data. Where the emission line is 
unresolved, the velocity width is an upper limit set by the effective width of the resolution element itself. 

f Red and blue side continuum measurements are variance— weighted averages made in 1200 A wide windows beginning 30 A from the 
wavelength of the peak pixel in the emission line. We employed a 10-iteration, 2a clipping algorithm to reduce the effect of spurious 
outliers occurring at long wavelength, where the sky noise is large. In some cases, a small correction factor was subtracted from the 
variance-weighted averages based on the detection of residual signal remaining in extractions of source-free, sky-subtracted regions of 
the two— dimensional spectra (see text, § 12.21 1. Quoted uncertainties account for photon counting errors in the source extractions added 
in quadrature to the photon counting errors derived in the blank-sky extractions. 

s 2<7 lower limit. The measurement of the red— side continuum for this source is formally consistent with no observable flux. The 
equivalent width limit was then set by using a 2a upper limit to f\ r in the expression given in footnote (c). 

h This line is unresolved. 



